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Abstract 

Using a recently developed method combining a nonspherical self-interaction corrected LDA+£7 
scheme and an on-site multi-body Hamiltonian [Phys. Rev. B 83, 085106 (2011)], we calculate the 
crystal field parameters and crystal field (CF) excitation levels of /-element dioxides in the fluorite 
structure with f n electronic configurations, including n — 1 (Pa02, PrC^), n = 2 (UO2), n = 3 
(Np02), and n — 4 (PUO2). It is shown that good agreement with experimental data (within 
approximately 10 to 20 meV) can be obtained in all cases. The properties of the multi-electron 
CF ground states are analyzed. 



PACS numbers: 71.70.Ch, 71.27.+a, 71.15.Mb 



I. INTRODUCTION 



The electronic structure of lanthanide and actinide compounds has a number of distinctive 
features that are manifestations of atomic /-electron physics in bulk solids, including strong 
on-site correlations and relativist ic spin-orbit effects. The effects of chemical environment 
on the ground states and excitation spectra of /-electrons are particularly interesting, since 
they are responsible for splitting the otherwise 2J+ 1-fold degenerate free-ion ground state 
2S+1 Lj^ giving rise to rich physics and applications. Recently, actinide dioxides in the cubic 
fluorite structure have attracted renewed theoretical interest in the context of their use as 
nuclear fueLP^. 

The crystal field (CF) method is a well established tool for describing the ligand environ- 
ment of localized electrons. In its conventional form, which requires spectroscopic informa- 
tion for fitting the CF parameters, the CF method has been applied to /-electron compounds 
with considerable succesd^Ei], including numerous characterizations of actinide oxides EH23I 
Using first-principles density functional theory (DFT) approaches, Divis and co-workers 
calculated the crystal field in praseodymium oxideiP^ID anc [ Colarieti-Tosti and co-workers 
studied PuO^S. However, since the CF splitting is much weaker that the Coulomb repulsion 
and spin-orbit coupling (SOC), DFT-based calculations are often plagued by various techni- 
cal issues, such as lack of a fully self-consistent treatment of the /-charge density or explicit 
consideration of electronic correlation. Recently Gaigalas and co-worker^l calculated CF 
levels of actinide dioxides with relativistic quantum chemical methods. 

We have recently developed a fully self-consistent method of calculating the CF parame- 
ters, which combines an improved nonspherical self-interaction free LDA+[7 scheme^ with 
a model on-site Hamiltonian including Coulomb, spin-orbit, and CF terms ™ Our approach 
utilizes the existence of multiple local minima in the LDA+U total energy functional and 
uses the corresponding /-electron wavefunctions and total energies to extract CF param- 
eters. Good agreement with experiment was obtained in terms of the predicted UO2 CF 
excitation spectrum (within about 10 to 20 meV) and magnetic properties of U02.^In this 
paper, we extend this method to calculate the CF parameters of other /-element dioxides 
MO2 in the fluorite structure with the f n configuration, including n — 1 (Pa02, Pr02), n = 3 
(Np02), and n = 4 (P11O2). Some results for UO2 (n = 2) are included for completeness. 
Other /-elements are not considered either because they have no valence / electrons (Ce02, 
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ThC^), no stable dioxides (heavier lanthanides), or no suitable pseudopotential presently 
available to us (Am02, CmC^). 



II. METHOD 

The CF of MO2 in the fluorite structure is given by: 

+ 32y^V 6 (V? - VllmcY*), (1) 

where V4 and Vq are CF parameters of the cubically coordinated metal ion, related to other 
common CF notations^ by 

V 4 = B4/8, 
V 4 = Be/16. 

In addition, free-ion parameters F k (k = 2, 4, 6) and ( describe the Coulomb and SOC terms, 
respectively, in the total Hamiltonian: 

H = H C F + Vee + (l'S- (2) 

where V ee designates the Coulomb repulsion between /-electrons. Since the Slater inte- 
grals F h (h = 2, 4, 6) in V^e are heavily correlatecpH, the following approximation has been 
adopted! 32 

F 2 = F 4 /0.668 = F 6 /0.494, (3) 

eliminating free parameters F 4 and F 6 . The Hamiltonian of Eq. ^ is diagonalized with f n 
basis wavef unctions, which are chosen in this work as n-body Slater determinants constructed 
from 14 f 1 spin-orbitals {Y™a} (m = — 3, •••,3, a =t, \). Therefore there are Gy± basis 
wavefunctions to expand an f n state. 

All DFT calculations were carried out using the same computational settings as in our 
previous workP^l. Input parameters for the LDA+lP^^ corrections are chosen as U = 6 eV, 
c = 0.5 and the J parameter for exchange interactions are determined by the requirement 
of numerical degeneracy of degenerate ionic states^. For each compound, 50 calculations 
with randomly initialized f n wavefunctions were carried out at the experimental lattice 
parameters (Table [I]). The magnetization axis for analyzing the energy eigenstates is chosen 
along z. More details of our technical approach can be found in Ref. |30l 
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TABLE I. List of studied /-element oxides M0 2 (M=Pr, Pa, U, Np, Pu), including the lattice 
constants a, number of localized /-electrons n, free-ion and crystal field ground states (GS) and 
their multiplicities (in parentheses), and the J parameter used in LDA+£7 calculations.^ 

III. RESULTS AND DISCUSSIONS 

The calculated model parameters are summarized in Table [TTJ For comparison, we also 
give the free-ion parameters F k and ( of tetravalent actinides in the corresponding fluorides.^ 
The Coulomb interactions F k do not enter the Hamiltonian of the f 1 compounds Pr02 and 
Pa02- F k is found to be slightly smaller in Np02 than UO2 and PUO2, in agreement with 
the trend observed in MF 4 . The SOC parameters (5/ « 0.2 to 0.3 eV of the heavier actinides 
are found substantially larger than lanthanide, since relativistic effects are more pronounced 
in heavier elements. The calculated (5/ of Pa02 is almost twice as large as the corresponding 
C4/ = 0.115 eV in the rare earth compound Pr02- C is predicted to increase over the actinide 
series, in agreement with experiment. However, our calculated (5/ values are overestimated 
by 5 to 15%. 

Higher localization of the 4/ states explains the smaller CF parameters V4 and Vq in Pr02 
compared to 5/ actinides. The 4th-order CF parameter V4 is significantly larger than Vq for 
all the dioxides, in agreement with results obtained from fitting experimental spectraP^. 

Using the parameters given in Table [TTJ the crystal field eigenstates are obtained by 
diagonalizing the effective Hamiltonian in Eq. [2| the resulting wave functions are visualized 
in Fig. [I] Following the procedure of Refs. [371 the radius R(Q = 9, cj)) of the spherical plots 
of the charge distribution is 

R{Sl) = (p(fi) - p) 1/3 , 
4 
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Free-ion parameters in MF4 from Ref. 1351 
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TABLE II. Calculated model parameters and energy eigenvalues in eV. F 4,6 may be derived from 
eq. |3j Excited CF levels are labeled with the corresponding degeneracy in parentheses and compared 
with reported measured or calculated values. 

where p(Q) is the spherical part of the charge distribution centered at the metal ion, and p 
is an appropriate amount of monopole subtracted from p(Q) to emphasize its asymmetric 
character. 

The predicted low-energy CF excitation levels are shown in Table [TT] alongside available 

5 



FIG. 1. The f n (n = 1, • • • , 4) eigenstates under the fluorite cubic crystal field. Only the 2 J + 1 
lowest states are shown. 



experimental data. Qualitatively, the correct ground states and ordering of the excited 
states are obtained in all cases. Quantitatively, good agreement with the measured spectrum 
has been obtained, with the errors in the excitation energies being within 10 to 20 meV. 
Reasonable agreement has also been found with previous theoretical calculations reported 
in Refs. [T9l [271 andl28l Next, we discuss each compound in detail (except UO2). 
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A. Pr0 2 and Pa0 2 



The f 1 compounds Pr02 and Pa02 differ from n > 1 cases in that the multiple local- 
minima issues that plague calculations for multi-/ electron systems are less severe. Out of 
the 50 random wave function initializations, approximately 5% with highly unstable starting 
states failed to converge within 100 electronic steps and were discarded. The rest exhibited 
a relatively small energy spread and were all within 0.2 to 0.3 eV from the CF ground state, 
compared to the spread of about 2 eV observed in U02-^ This shows that the many-body 
interaction is a main reason for the existence of many local-minimum solutions, and without 
this obstacle the f 1 calculations can find the j = 5/2 Russell-Sanders ground state, even 
though they may fail in locating the CF ground state. The original LDA+[/ 34 scheme was 
tested for Pr02 and Pa02 and found to increases significantly the energy spread of the 
local-minimum solutions due to orbital-dependent self-interaction errors 29 . 

The predicted T 8 —> T 7 excitation energy for Pr02 is 129 meV, in excellent agreement with 
the measured value of 131 meV^ and more accurate than our previous rough estimation of 
73 to 142 meV in Ref . [29l showing that our method based on Eq. ^ leads to significant error 
cancellation in the calculated CF energies. Predictions for the higher CF levels of J = 7/2 
are 0.376 (r' 6 ), 0.433 (T f 8 ) and 0.622 eV (r' 7 ), respectively, compared with observed values of 
0.320, 0.390 and 0.580 eV from Ref. [23l The relative splitting within the J = 7/2 manifold 
agrees very well with experiment, showing the validity of our predicted CF parameters, while 
the center of these levels are 11% too high, due to the over-estimated spin-orbit coupling 
(our ( = 0.115 eV compared to 0.1 eV of Ref. 1231) . Experiments on Pa02 are relatively 
scarce. The only available number of 140 meV for the r 8 T 7 transition cited in Ref. |38] 
is based on private communications, which we inquired about but could not confirm. Our 
prediction of 186 meV for T 8 —> T 7 in Pa02 is substantially larger than the corresponding 
value for Pr02, in agreement with the trends in CF parameters in Table |TT} 



B. U0 2 

Recent measurement by Nakotte et a/.^of the crystal field levels in UO2 provides updated 
information than Amoretti et a/. 15 : the excitation peak at 180 meV is spurious. 
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(5) r 7 ® 
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0.10 
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(6) r 7 ® 


-0.06 




-0.73 


-0.10 




-0.70 



TABLE III. The Tg ground state quartet and I\ excited doublet of Pr02 and Pa02 (f 1 ) in the 
fluorite structure, and the corresponding spin and total magnetic moment in jig. 



C. Np0 2 

Np0 2 has the 5/ 3 configuration. The excitation energy between the CF ground state 
and the first excited state has been measured to be 55 meV 16 . Our prediction of 34 
meV is a reasonable under-estimation. Note that a recent quantum chemical calculatioiP^ 
for Np0 2 predicted the same value as ours. Two estimated values for the second excited T 6 
energy level (145 meV and 274 meV) are given in Ref. [16j and only the first value scales over 
the actinide dioxides series 19 . We predict an excitation energy of 125 meV for Tg 2 ^ — >> T 6 , 
which agrees with the latter assessment. 

Since the f 1 configuration in the fluorite structure is split by the crystal field into the T 8 



quartet and T 7 doublet (see Table III), the multi-electron configurations f n of U0 2 , Np0 2 
and Pu0 2 are sometimes interpreted within a picture where the added electrons gradually 
fill the CF levels, in analogy to the well-known scenario of d-electrons filling the £ 2 ^ and 



e g CF levels in transition metal compounds. As we have shown previously 30 , this picture 
fortuitously holds for the T 5 ground state of U0 2 . However, transition metal CF splittings 
are usually several eV, while for / the CF splittings (on the order of 0.1 eV) is much smaller 
than the effective Coulomb interactions (~ eV). Hence, the f n eigenstates are in general 



(2) 

multi-configurational. According to Table IV, the Tg ground states of Np0 2 are composed 



of multiple determinants, including ones with substantial projections onto not only the T 8 
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A* 


Proj. 


«€ 


-1.00 


1.72 


0.68(1,2,3) + 0.57(1,5,6) 




-0.37 


0.57 


0.68(1,2,4) + 0.57(2,5,6) 




0.37 


-0.57 


0.68(1,3,4) + 0.57(3,5,6) 




1.00 


-1.72 


0.68(2,3,4) + 0.57(4,5,6) 



TABLE IV. The Ig ' ground state quartet of Np02 (5/ 3 ), along with their magnetic moment and 
projection onto Slater determinants (only a few leading terms shown) composed of f 1 eigenstates 
Ts and IV One-electron orbitals (1) through (6) are defined in Table [III} 

CF ground states, but also the T 7 excited states of f 1 . In other words, the f 3 ground state 
Ig occupies both T 8 and T 7 orbitals in order to lower its electrostatic energy at the expense 
of a slightly increased CF energy. 



D. Pu0 2 

Pu02 has the 5/ 4 configuration. The crystal field was measured by Kern and co-workers 
using inelastic neutron scattering (INS) 17 . Our calculated Ti —> T 4 excitation energy of 97 
meV agrees reasonably well with the measured value of 123 meV 17 and a previous calculation 
of 99 meV by Colarieti-Tosti et al 27 . Note that the splitting was under-estimated in all the 
calculations, including this work and Refs. HU [271 an d EHl The non-magnetic 5/ 4 ground 
state Ti with ji$ — \i — 0, 

& = 0.7(1, 2, 3, 4) + 0.32(1, 4, 5, 6) + 0.32(2, 3, 5, 6) + ... , 

is sometimes referred to as four fully filled T 8 f 1 orbitals. We obtained | (Ti 1 1 , 2, 3, 4)| 2 = 
0.49, showing that such a simplified picture of 4 filled T 8 orbitals is not entirely valid and 
multi-electron correlations account for more than 50% of the ground state wave function. 

As a simple application of the CF calculations, Fig. [2] shows the calculated electronic 
entropy 

£ e = -J^ Pi In pi, 

i 

where pi = e -^A# T / ^\ e -Ej/k B T - g ^ e Boltzmann probability of the electronic eigenstate 
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i. As shown by Konings, 38 the vibrational contribution to the total entropy of actinide 
oxides varies smoothly across the elemental series, while electronic contributions, which 
depend delicately on the CF excitation energies, cannot be interpolated over the series. 
To accurately predict thermodynamic properties of actinide oxides, the electronic entropy 
cannot be ignored. Our predicted S e (solid curves) agree reasonably well with the results of 
Ref. 38 (crosses) at T = 298.15 K. 

In conclusion, we have calculated the CF levels of Pr02, Pa02, Np02, and Pu02. The f- 
electron charge density and on-site correlations are calculated fully self-consistently within a 
version of LDA+U that removes orbital-dependent self-interaction energies. Good agreement 
with experimental CF levels and a consistent trend across the actinide series have been 
achieved. In both Np02 and Pu02, substantial contributions of the T 7 one-electron excited 
state are found in the multi-electron crystal field ground states. 
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